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Abstract 

Heterodyne laser phase measurements in a space-based gravitational wave interferometer are 
degraded by the phase fluctuations of the onboard clocks, resulting in unacceptable sensitivity 
performance levels of the interferometric data. In order to calibrate out the clock phase noises 
it has been previously suggested that additional inter-spacecraft phase measurements must be 
performed by modulating the laser beams. This technique, however, considerably increases system 
complexity and probability of subsystem failure. 

With the advent of self-referenced optical frequency combs, it is possible to generate the hetero¬ 
dyne microwave signal that is coherently referenced to the onboard laser. We show in this case that 
the microwave noise can be cancelled directly by applying modified second-generation Time-Delay 
Interferometric combinations to the heterodyne phase measurements. This approach avoids use of 
modulated laser beams as well as the need of additional ultra-stable oscillator clocks. 
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I. INTRODUCTION 


Gravitational waves (GWs) are predicted by Einstein’s theory of general relativity and 
represent disturbances of space-time propagating at the speed of light. Because of their 
extremely small amplitudes and interaction cross-sections, GWs carry information about 
regions of the Universe that would be otherwise unobtainable through the electromagnetic 
spectrum. Once detected, GWs will allow us to open a new observational window to the 


Q. 


Universe, and perform a unique test of general relativity 

Since the hrst pioneering experiments by Joseph Weber in the early sixties j^, several 
experimental groups around the world have been attempting to detect GWs. The coming 
on-line of next-generation gravitational wave interferometers j^, however, is likely to 
break this sequence of experimental drawbacks and result into the hrst detection before the 
end of the current decade. 

Gontrary to ground-based detectors, which are sensitive to gravitational waves in a band 
from about a few tens of Hz to a few kilohertz, space-based interferometers are expected to 
access the frequency region from a few tenths of millihertz to about a few tens of Hz, where 
GW signals are expected to be larger in number and characterized by larger amplitudes. 
The most notable example of a space interferometer, which for several decades has been 
jointly studied in Europe and in the United States of America, is the Laser Interferometer 


Space Antenna (LISA) mission j^. By relying on coherent laser beams exchanged among 
three remote spacecraft along the arms of their forming giant (almost) equilateral triangle of 
arm-length equal to 5 x 10® km, LISA aims to detect and study cosmic gravitational waves 
in the 10“^ — 1 Hz band. 

A space-based laser interferometer detector of gravitational waves measures relative fre¬ 
quency changes experienced by coherent laser beams exchanged by three pairs of spacecraft. 
As the laser beams are received, they are made to interfere with the outgoing laser light. 
Since the received and receiving frequencies of the laser beams can be different by tens 
to perhaps hundreds of MHz (consequence of the Doppler effect from the relative inter¬ 
spacecraft velocities and the intrinsic frequency differences of the lasers), to remove these 
large beat-notes present in the heterodyne measurements one relies on the use of a mi¬ 
crowave signal generated by an onboard clock (usually referred to as Ultra-Stable Oscillator 
(USO)). The magnitude of the frequency fluctuations introduced by the USO into the het- 
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erodyne measurements depends linearly on the USOs’ noises themselves and the heterodyne 
beat-note frequencies determined by the inter-spacecraft relative velocities. Space-qualihed, 
state of the art clocks are oven-stabilized crystals characterized by an Allan deviation of 




10 


-13 


for averaging times of 1 — 


terest to space-based interferometers 


000 s, covering most of the frequency band of in- 
7|. In the case of the LISA mission, in particular, 


it was estimated js] that the magnitude of the power spectral density of the USO’s relative 
frequency fluctuations appearing, for instance, in the unequal-arm Michelson Time-Delay 
Interferometry (TDI) combination X would be about six orders of magnitude larger than 
those due to the residual (optical path and proof-mass) noise sources. 

A technique for removing the USO noise from the TDI combinations was devised (see 


Q 


for more details). This technique requires the modulation of the laser beams exchanged 
by the spacecraft, and the further measurement of six more inter-spacecraft relative phases 
by comparing the sidebands of the received beam against sidebands of the transmitted beam. 
The physical reason behind the use of modulated beams for calibrating the USOs noises is 
to exchange the USOs phase fluctuations with the same time delay as their lasers among 


the three spacecraft by performing side-bands/side-bands measurements [SHlOj]. In so doing. 


additional six phase measurements are generated that allow one to calibrate out the USOs 
phase fluctuations from the TDI combinations while preserving the gravitational wave signal 
in the resulting USO-calibrated TDI data. 

It should be noticed, however, that if we could coherently transfer the laser phase fluctu¬ 
ations to the microwave signal used in the heterodyne measurements, then we would need to 
cancel only one noise (the laser noise), which might be possible by deriving some new TDI 
combinations. Coherently linking optical laser frequencies to microwave frequencies has been 
thought to be impossible because of the inability to directly count the optical frequency of a 
laser. However, with the recent advent of the self-referenced octave-span optical frequency 
comb (OFC) scheme (for which Hall Hj] and Hansch 12(] received the physics Nobel Prize 
in 2005) it is now possible to generate a microwave frequency signal that is coherent to the 
frequency of the laser at a level signihcantly better than the frequency stability required of 
a USO to avoid the modulation-driven USO noise calibration procedure. 

Although most of the reported implementations of the self-referenced frequency comb 
technique have been performed in a labora tory environment, active developments are being 
made for space-qualihed OFC systems 13|, [l^ . In addition, more recent micro comb source 
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developments promise much smaller and integratable comb devices 1^. If each spacecraft 
could then generate heterodyne measurements by relying on the OFC technique, then the 
question to be answered is whether there exist new TDI combinations that could account for 
the modihed laser noise transfer functions into the heterodyne measurements. The answer 
to this question is yes, and in what follows we will derive their expressions. 

In Section mi we give a brief summary of the OFC technique, and of characteristics of 
a self-referenced comb. In Section m we then exemplify through a simple interferometric 
conhguration illustrating how the “USO noise” can be canceled with TDI when the OFC 
technique is implemented. In the same section we then move on to the realistic conhguration 
of a space-based interferometer (such as LISA) and derive the expressions for the twelve 
inter-spacecraft heterodyne phase measurements when the OFC technique is used. This 
allows us to obtain the new second-generation (i.e. hex-free) TDI combinations 16l-ll9| that 
simultaneously cancel the laser and the comb-generated microwave signal phase noises. In 
Section IIVI we hnally present our concluding remarks, and emphasize that the use of the 
OFC technique will result into a hardware and system design simplihcation, and increase 
system reliability of future space-based gravitational wave interferometers. 


II. THE OPTICAL FREQUENCY COMB TECHNIQUE 

An OFC is a set of narrow spectral lines equally spaced in the frequency domain. More 
importantly, when all the comb lines are phase coherent (i.e. have hxed phase relationship) 
the superposition of all the comb lines manifest themselves as a periodical wave train in the 
time domain. The resulting repetition rate is equal to the comb spacing and it is governed 
by the Fourier transform relationship. More often, the comb lines are in phase so the time 
domain waveform is a train of short pulses separated by the inverse of the repetition rate. 
The pulse width is then Fourier transform limited by the comb spectral width. 

Typical OFCs are generated by mode-locked lasers in which multiple modes are locked in 
phase, outputting short pulse trains. The pulse train propagates with a group velocity while 
the underline optical carrier travels at the phase velocity. When the modes are exactly in 
phase, the peak of the carrier is aligned with the peak of the pulse envelope, as shown in the 
hrst pulse in hgure [TJ In the absence of any dispersion, this relationship holds over time 
as the pulse propagates. In practice, dispersion changes the relative phases of the modes. 
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and the peak of the carrier amplitude shifts from that of the pulse envelope. This offset is 
typically referred to as Carrier Envelope Offset (CEO) phase shift and the rate of this offset 
phase change over a period gives the CEO frequency 2^, as shown in the frequency domain 
spectrum in hgure [H In particular, this implies the well known mode frequency relation: 
fm = fceo + T^frep- Siuce fceo IS affected by dispersion, as there are many environmental 
effects influencing it, fceo is not a stable parameter and, as a consequence, there is no hxed 
frequency or phase relationship between the optical carrier and the pulse envelope, i.e. the 
beat-note of the rf signals among the modes. For many years before the recent self-referenced 
OFC development, the two frequencies were kept separate and used independently. 

Recent advancements in the study of ultra-fast phenomena and the held of optical fre¬ 
quency metrology have revolutionized the use of OFC with abilities to measure and control 
the CEO. As shown in reference 20j, if one compares the A^mth harmonic of the mth mode 


with the A"rn'th harmonic of the m'th mode, the beat note of the two yields {Nm — Nm')fceo- 
In particular, if the comb has an octave spanning width where one can choose = 2 and 
Nm' = 1, the mixing beat-note of lf-2f directly gives the measure of the CEO frequency. 
Like in a typical phase locking loop, the measured CEO phase can be used to stabilize the 
CEO itself, resulting in an OFC that has the well-dehned relationship between the optical 
frequencies fm in the modes and the rf frequency of the repetition rate frep, now referred to 
as self-referencing j^|. Such control capability not only allows one to align carrier amplitude 
peaks with the short pulse envelopes in ultrafast phenomena study, but more revolutionarily 
provides the optical frequency divider ability down to microwave in one simple step. In the 
latter case with a self-referenced comb, what it means is that the optical held in the comb 
has a coherent phase relationship with its rf beat-note signal and one can count the optical 
frequency by counting the rf frequency. Furthermore, the optical phase noise is also directly 
down converted to the rf phase noise. Recent studies have shown that the comb frequency 
precision can be at 1 x 10“^® 22| while the CEO phase can be controlled down to the 


mrad level. In the context of the investigation in this paper, the properties of the phase 
coherent optical divider are exploited for the onboard clock phase noise cancellation and for 
the generation of stable onboard rf signals directly from the comb. 

The most mature self-referenced OFCs are generated with a combination of femtosecond 
mode locked laser and super-continuum generation to achieve octave spanning spectrum 
width. While both Ti:Saphire and hber mode locked lasers are successfully implemented for 
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FIG. 1: Illustration of a mode locked laser pulse train and its corresponding comb in the frequency 
domain. The top plot in the time domain shows the carrier envelope offset phase and its changes. 
The bottom plot shows the comb lines and the carrier envelope offset frequency near the DC 
frequency. 

self referencing, the most snitable laser sonrce for space interferometer mission may be fiber 
laser based systems. In fact, self-referencing comb systems are now commercially available, 
mostly for laboratory applications. Efforts are being made in developing the fiber OFC 
technology for space applications. A flight OFC system is being developed for sonnding 
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rocket experiments [13[ , and a femtosecond fiber laser has already been flown in space . 
It is in fact expected that a flight OFC system, suitable for space missions, will become 
available in the near future. 

More recently it has been shown that coherent OFC can also be generated using high-Q 
micro resonators through Kerr nonlinear four wave mixing process 23|. In this case, a cw 
laser is resonantly injected into the resonator and strongly conhned in the small volume 
mode. The high Q-factor leads further held buildup in the resonator, which enhances the 
Kerr nonlinear effect in the micro resonator. Various cavity eigenmodes can, therefore, be 
excited through four-wave mixing as the laser light sequentially cascades from the pump to 
the other modes 2^, forming an optical comb with spacing determined by the resonator 


mode space, free spectral range. When the condition is right, these modes are phase locked 
in a soliton form, very much like a typical mode locked laser 2^. In comparison with their 


mode-locked lasers counterparts, however, these whispering galery mode optical frequency 
comb generators are characterized by a signihcantly reduced size and power consumption, 
along with a high repetition rate. They are, therefore, particularly suitable for miniaturiza¬ 
tion, chip integration, and space applications. For the purpose of self-referencing in optical 
metrology, it is important that the comb spans over one octave for lf-2f self-reference, or 
at least 2/3 of an octave for 2f-3f self-referencing. While self-referenced micro comb are 
still been developed in research labs, we believe space-qualihed, compact Kerr combs will be 
available in the near future. 

The striking property of the self-referenced OFC is the phase coherence between the 
optical carriers and their rf beat note signal. As it will be shown later, this allows the laser 
and resulting microwave signal phase noises cancellation elegantly in the generalized TDI 
combinations. Furthermore, one can take advantage of the onboard frequency stabilized 
laser system to generate a stable microwave signal. Since the onboard stabilized lasers are 
expected to achieve a stability at the level of 10 “^“^/\/Hz or better over the frequency range 
of interests to space interferometer gravitational wave detection, by phase-locking the optical 
frequency comb to the stabilized laser, the corresponding rf beat-note stability will be also 
of a few parts in 10^^. Therefore, the stabilized OFC technique will provide an rf signal 
that can be signihcantly more stable of the current best space-qualihed USOs. The OFC 
implementation will not only eliminate the need for separate USOs but also remove entirely 
the need for modulated beams used by the USO calibration method, as it will be shown in 
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the following sections. 


III. TIME-DELAY INTERFEROMETRY WITH OPTICAL FREQUENCY COMB 


In the following subsections we derive the modihed expressions of the second- 
generation ly, ll8|, ll9| TDI combinations that simultaneously cancel the laser and the mi¬ 
crowave phase noises when the OFC technique is implemented. To physically understand 
how this is possible, we will hrst consider a simple interferometer conhguration in which the 
rate-of-changes of the arm-lengths are large enough to require the removal of the beat-notes 
from the two-way data but small enough to treat the arm-lengths as constant. This will 
allow us to establish the basis for then deriving new second-generation TDI combinations 
that a three-arms, six laser links space-based interferometer can generate. 


A. TDI with OFC: an example application 

Let us consider the following idealized unequal-arm interferometer whose end-mirrors are 
moving with velocities Vi and V 2 along the directions of arms 1 and 2 respectively (see 
£gure[2]). These velocities are relative to the inertial reference frame in which the laser and 
the beam-splitters are at rest. The two light beams coming out of the two arms are not 
made to interfere at a common photodetector as each is made to interfere with the incoming 
light from the laser at a photodetector. This is because direct recombination of the beams 
at a single photodetector would not suppress the laser phase noise in the phase difference 
measurement due to the inequality of the arms. ^ 

The velocities of the mirrors can be assumed to have frequency components that are 
outside the band of operation of the interferometer, while their amplitudes are such as to 
require the use of a numerically controlled oscillator (NCO) (driven by a clock) to remove 
the beat-notes in the two “two-way” heterodyne phase measurements, 0i(t) and 02 (i)- As a 

^ By using two independent photodetectors we decouple the phase fluctuations experienced by the two beams 
in the two arms and obtain two (rather than one) “two-way” measurements. This “data redundancy” 
allows us to cancel the laser phase noise in the data while retaining a possibly present gravitational wave 
signal. This is done by properly time-shifting and linearly combining in digital form the two two-way 
phase differences. 


8 









Vo,P (t) 


yi(t) 


FIG. 2: Light from a laser is split into two beams, each injected into an arm formed by pairs 
of mirrors moving with velocities Vi and V 2 relative to the inertial reference frame with respect 
to which the laser and the beam-splitter are at rest. Since the length of the two arms, Li and 
L 2 , are different, the light beams from the two arms are not recombined at one photo detector. 
Instead each is separately made to interfere with the light that is injected into the arms. Two 
distinct photo detectors are now used, and phase (or frequency) fluctuations are then monitored 
and recorded there. See text for details. 

result of this “down-conversion” operation, the two-way data assume the following forms ^ 

0i(^) = - 2Vi) - 1^0 - aif]t + p(t - 2Li)-p(t) - aig(t) + hi(t)+ni(t), (1) 

02(t) = K(1 - 2 V 2 ) - z/Q - a2f]t + pit - 2 L 2 ) - pit) - a2qit) + / 12 (f) + n2(t), (2) 

where ai, 02 are selected by the NCO in such a way to suppress the beat-notes from the 

^ Throughout this article the speed of light c and 27r are assumed to be equal to 1 
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two phase measurements to a required value, / is the microwave frequency generated by 
the clock driving the down-conversion, hi, h 2 are the contributions from a possibly present 
gravitational wave signal, and ni, n 2 are the random processes associated to the overall 
remaining noises affecting the two phase measurements. Note that in Eqs. ([H |2]) we have 
denoted with q{t) the random process associated with the phase noise from the clock driving 
the NCO |8| and p{t) the laser noise. Since the coefficients a* , i = 1, 2 can be made to be 
equal to = —2ViZ/o// , i = 1,2, then Eqs. ([H [2]) become 

4>i{t) = p{t - 2Li) - p{t) - aiq{t) + hi{t) + ni{t) , (3) 

02(t) = pit - 2 L 2 ) - pit) - a2qit) + h2it) + n2it) , (4) 

If we now assume the clock frequency / to be the result of implementing a self-referenced 

OFC that is driven by the onboard laser, then the following relationship between the phase 

fluctuations of the local oscillator, g(f), and those of the laser, pit), holds 


q[t) = ^ pit) + Aqit) . 
^0 


(5) 


In Eq. (IT^ we have denoted with Ag(f) the residual noise representing the level of coherence 
between the frequency of the laser and the microwave frequency referenced to it. It has been 
demonstrated experimentally 22| that the magnitude of this residual phase noise can be 
reduced to a level corresponding to an Allan standard deviation of a few parts in 10“^® over 
time scales of interest to space-based GW interferometers. Since this value is three orders of 


magnitude smaller than the USO noise 
procedure based on modulated beams 


evel required for avoiding the use of the calibration 
Sj, we can then safely disregard it in our analysis. 
After substituting Eq. ([5]) into Eqs. (El 0]) we obtain the expressions for the two-way phase 
differences when the OFC technique is implemented 


01 (f) = pit - 2Li) - pit) + 2Vipit) + hiit) + niit) , (6) 

02(t) = pit - 2 L 2 ) - pit) + 2V2pit) + h2it) n2it) . (7) 


Following 


26| it is then straightforward to show that the following combination of the two 


two-way phase measurements cancels the laser noise 


X^^^it) = [01 (f - 2 L 2 ) - (1 - 2^2)01 (f)] - [hit - 2Ei) - (1 - 2Vi)hit)] . (8) 


In other words, all pit) and g(f) terms (through Eq. [5]) drop out in . Although the 

above expression of the unequal-arm Michelson TDI combination reflects the assumption 
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that the velocities of the mirrors are large enough to require the removal of the beat-notes 
from the two-way data but small enough to treat the arm-lengths as constant, it shows 
that the implementation of the OFC technique allows us to simultaneously cancel the laser 
and microwave noises with TDI. In the following sections we will remove the assumption of 
constant arm-lengths and derive the TDI expressions that account for the rotation of the 
array (Sagnac effect) as well as the time-dependence (“flexi”) of its arms 


16 


18 


19|. 


B. Generalized TDI formulation 


In order to derive the new TDI combinations valid when the OFC technique is used for 
generating a microwave signal coherent to the onboard laser, we will follow the description 
of a space-based interferometer (such as LISA and eLISA/NGO with its three operational 
arms js-?]) adopted in earlier publications (see reference 0 and references therein). 

As shown in Figure |3l where the overall geometry of the interferometer is dehned, there 
are three spacecraft, six optical benches, six lasers, six proof-masses, and twelve photo 
detectors. There are also six phase difference data going clock-wise and counter-clock-wise 
around the triangle. The spacecraft are labeled 1,2,3 and the optical paths traveled by 
the light beams as they propagate along the arms of the constellation are denoted Lj and 
Lj/, with z = 1, 2, 3 being opposite spacecraft i. The two optical paths along the same arm 
are different because the rotational motion of the array results into a difference of the light 
travel times in the two directions around a Sagnac circuit jl6, 27|. As shown in Figure [3] we 
have conventionally denoted with Lj/ and Li the optical delays for clockwise and counter- 
clock-wise propagation respectively. Furthermore, since Li and Lj/ not only differ from one 
another but can be time dependent (they “flex”), new second-generation TDI combinations 
were derived in order to account for this effect. This was done by using the non-commuting 
time-delay operators formalism discussed in [^, and we will rely on it for the derivation of 
the new second-generation TDI combinations that simultaneously cancel the noises of the 
microwave signals generated by the OFC technique. 

In Figure [3] the vertices 1,2,3 are oriented clockwise, while the unit vectors between 
spacecraft are denoted hj, oriented as indicated in the hgure. We index the phase difference 
data to be analyzed as follows: the beam arriving at spacecraft i has subscript i and is primed 
or unprimed depending on whether the beam is traveling clockwise or counter-clockwise (the 
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FIG. 3: Schematic array configuration. The spacecraft are labeled 1, 2, and 3. The optical paths 
are denoted by Li, Li/ where the index i corresponds to the opposite spacecraft. The unit vectors 
hj point between pairs of spacecraft, with the orientation indicated. 


sense defined here with reference to Fignre[3]) aronnd the triangle, respectively. Thns, as seen 
from the hgnre, Si is the phase difference time series measnred at reception at spacecraft 1 
with transmission from spacecraft 2 (along L 3 ). 

Similarly, si/ is the phase difference series derived from reception at spacecraft 1 with 
transmission from spacecraft 3. The other fonr one-way phase difference time series from 
signals exchanged between the spacecraft are obtained by cyclic permntation of the indices: 
1 ^ 2 —)■ 3 —)■ 1. We also adopt the notation for delayed data streams described in [l9|, in 
which there are six time-delay operators "Dj , Pj/, i = 1,2, 3, and where, for any data stream 

- Li) , Vi,^{t) = - Li>) . (9) 


For example, 'D 2 Si{t) = si(t —L 2 ), 'D 2 'T>ySi{t) = si{t — Ly {t — L 2 ') — L 2 '(t)), etc. Note that 
the operators do not commute, as: VyVySiit) = si(t — Ly (t — Ly) — Ly (t)) 7 ^ VyVySiit) = 

Si{t — Ly{t — L 2 '{t)). 

Following jsl, the heterodyne phase measurements, Si{t) , s*/ , i = 1,2,3 can be written 
in the following form 
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Si{t) = [i^2'(l - Ls) -Vi- aifi]t + V^P 2 ' - Pi{t) - aiqi{t) + + Ni{t) , 

sy{t) = - L 21 ) - vv - ai>fv]t+ V 2 ip^-Pv{t) - avqv{t) + + Ny{t) , (10) 


where we have denoted with pi, pi' the lasers phase fluctuations, the microwave frequencies 
/i , fy are generated by the OFCs driven by lasers 1 and 1' respectively, and the coefficients 
Oi, ay are syntl^sized by a numerically-controlled oscillator (NCO) to be equal to the 
following values 


Z/2/(l - L3) - Z/i U^{l-L2')-Vy 

( 1 \ — - - - , (l\f — - 


( 11 ) 


fl ’ fv 

In Eq. (fTOj) we have denoted with Ni , Ny all the remaining noises affecting the one-way 
measurements, and with , sf^ the contribution to the one-way Doppler measurements 
from a passing gravitational wave signal. 

Consistently with earlier work |8|, Eqs. (1T0|1 reflect the assumption of having the fre¬ 
quencies of the lasers different from each other. Note that the phase fluctuations from the 
microwave frequency references driving the NCO enter into Si and Si/ through the terms a^qi 
and ayqy respectively. Since these phase fluctuations associated to the microwave signals 
have been generated by the frequency comb subsystems, they are related to the laser phase 
fluctuations, pi and py, through the following relationships 


gi(t) = —Pi(i) + Agi(f) , qy(t) = —py(t) + Aqy(t) . 

Ui Uy 


( 12 ) 


In Eq. flT^ we have denoted with Aqi{t), Aqy{t) the residual noises representing the level 
of coherence between the frequencies of the lasers and the microwave frequencies referenced 
to them. The magnitude of this residual phase noise has been demonstrated experimentally 


22l | to correspond to an Allan standard deviation of a few parts in 10 over time scales of 


interest to space-based GW interferometers. Since this is three orders of magnitude smaller 
than the USO noise required for avoiding the use of the calibration procedure based on 
modulated beams j^, we infer that the contribution from the noises Ag* , Aqy , i = 1, 2, 3 
to the overall noise budget of the TDI combinations can be regarded as negligible. For this 
reason from now on they will be disregarded in our analysis. 

If we now substitute Eqs. lfTTl [12]) into Eq. dTOj) . the heterodyne measurements Si, Si/ 
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FIG. 4: Schematic diagram of proof-masses-plus-optical-benches onboard the spacecraft. The left- 
hand bench reads out the phase signals si and ri. The right-hand bench analogously reads out si' 
and T\i. The random displacements of the two proof masses and two optical benches are indicated 
(lower case dj, (5j/ for the proof masses, upper case Aj, Aj/ for the optical benches). 


assume the following forms 

■Si(^) = - Ti(^) - AiPi{t) + sf^{t) + Ni{t) , 

sr{t) = T>2'P3-pv{t) - Ai>pr{t) -f + Ni>{t) , (13) 


where the coefficient Ai , Ay are equal to 
Oi/i 


Ai = 




z/2t(l - L3) - z/i _arfr 

- , Ay = 




Uy 


Z/3(l - La') - Uy 
Uy 


(14) 


Note that, over periods of a few light-times, the coefficients Ai, Ay {i = 1, 2, 3) can be 
regarded as constant since the spacecraft relative velocities start to display a noticeable 
change only over periods of weeks to months j^ |. 

If we assume the optical bench design to be as shown in Figure IH then six more phase 
difference series result from laser beams exchanged between adjacent optical benches within 
each spacecraft; these are similarly indexed as r*, r,/, i = 1,2,3. ^ The photo receivers that 
generate the data si, Si/, ri, and Ty at spacecraft 1 are also shown in hgurelH 


^ The optical bench design shown in figure |4] represents one of the possible configurations for integrating 
the onboard drag-free system with the heterodyne measurements. Although other optical bench designs 
result into different inter-proof-mass phase measurements [d, 29|, the resulting expressions for the new 
TDI combinations derived in this article are all equivalent. 
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The adopted optical bench design shows optical fibers transmitting signals both ways 
between adjacent benches. We ignore time-delay effects for these signals and will simply 
denote by fii{t) the phase fluctuations upon transmission through the fibers of the laser 
beams with frequencies z/*, and z/j/. The phase shifts within a given spacecraft might not 
be the same for large frequency differences z/, — z/j/. For the envisioned frequency differences 
(a few hundred MHz), however, the remaining fluctuations due to the optical fiber can be 
assumed to be independent of the direction of light propagation through them 30[ |. 

If we limit our attention only to the terms containing the laser phase noises, we can then 
write down the following expressions for the Tj, Tj/ phase measurements 


n{t) = + (15) 

f ^ 

Ty{t) = [z/i - Z/i/ - Ci>fi>]t + pi{t) -Pl'{t) - Cy—Pl'{t) , (16) 

Z/i/ 


where, like in the case of the inter-spacecraft measurements Si, Si/, the coefficient Ci, Ci/ are 
determined by the NCO to be equal to 


Cl = 


Z/i/ - z/i 


Cl' = 


Z/i - Z/i' 


(17) 


/i /i' 

Note that in Eq. flTTl) we have used the relationship between the phase noise of the microwave 


signal and that of its driving laser (Eq. ([12])). As discussed in other TDI references [sj, [19 1 


the inter-bench phase measurements r*, Tj/ enter into the TDI combinations in the form of 
their differences. It is therefore convenient to rewrite the set of heterodyne measurements 
in the following form (for simplicity of notation we will omit the contribution from the 
gravitational wave signal and the other noises in the si and si' measurements) 


Si 

Sl' 

^1 


T>3P2' - {1 + Ai)pi , 


T^2'P3 ~ (1 + Aii)pi' , 

D-D' /, , pv. 

= (1 + Y)ft' 


(i + y)pi. 


where pi, pi', are equal to 


Z/i' - z/i 



Pi' 


Z/i - Z/i' 


Z/i' 


(18) 

(19) 

( 20 ) 


( 21 ) 


Six other relations, for the readouts at vertices 2 and 3, are given by cyclic permutation of 
the indices in Eqs. flTSl fl^ l20|) . 
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In order to derive the new TDI combinations that rely on the OFC techniqne, we will 
start by taking specihc combinations of the one-way data entering in each of the expressions 


derived above. Following the approach in 1^, [lOj], these combinations are identihed by 


requiring the laser noises entering into them to be only of one “kind”, i.e. with primed or 
unprimed indeces. After some simple algebra to account for the presence of the heterodyne 
coefficients, these combinations of the one-way heterodyne measurements (called rji , rji/ in 
here assume the following forms 


1 + ^ ^ 3^2 _ 

Vl = , , I Si - = V3P2 - hpi , 


1 + £2 


1 _L £2 
^2 


1 + Ay 

Vr = si/ -h . = Pa'Ps - h'Pi , 

where the coefficients Ji , ly are given by the following expressions 

(l + AiO(l + f) n^zl.^n^£2: 


( 22 ) 

(23) 


h = 


1 + ^ 
^2 


_(1 + Ai)(l + ^) 

ly = — 


1 _L £2 
^2 


(24) 


and the other r^-combinations are obtained as usual by permutations of the spacecraft inde¬ 
ces. 

The new TDI combinations are chosen in such a way so as to retain only one of the three 
noises pi, i = 1,2,3, if possible. In this way we can then implement an iterative procedure 


based on the use of these basic combinations and of time-delay operators, to cancel 


noises after dropping terms that are quadratic in L/c or linear in the accelerations 
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le laser 


18 


19|. 


This iterative time-delay method, to hrst order in the velocity, is illustrated abstractly as 


follows. Given a function o 
standard comma notation 


time T = 4/(t), time delay by L, is now denoted either with the 


3l[ | or by applying the delay operator "Dj, 


D,T = T,i = T(f-L,(t)). 


(26) 


We then impose a second time delay Lj{t): 


VpA> = - Lj{t) - Li{t - Lj{t))) 

~ T(t - Lj{t) - Li{t) Li{t)Lj] 


(26) 
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A third time delay Lk(t) gives 


VkVjVi^ = - Lk{t) - Lj{t - Lk{t)) - Li{t - Lk{t) - Lj{t - Lk{t)))) 




,ijk 




,ijk 


Li{Lj + Lfc) + LjLk 


(27) 


and so on, recnrsively; each delay generates a first-order correction proportional to its rate 
of change times the snm of all delays coming after it in the snbscripts. Commas have now 
been replaced with semicolons to remind ns that we consider moving arrays. When 
the snm of these corrections to the terms of a data combination vanishes, the combination 
is called flex-free. Finally, it shonld be noticed that each delay operator, Pj has a nniqne 
inverse operator, whose expression can be derived by reqniring that V~^Vi = / and 

neglecting qnadratic and higher-order velocity terms. Its action to a time series '^{t) is 


Vr^^{t) = ^{t + L,{t + L,)) , 


(28) 


i.e. it advances the time-series by a delay Lj estimated not at time t but rather at time 
t Li. 


C. The unequal-arm Michelson 

The unequal-arm Michelson combination synthesized onboard spacecraft ^ 1 relies on 
the four measurements rji, rjii, r] 2 ', and From Eqs. fl2^ [23]) . after some simple algebra, it 
is easy to show that the following two combinations, 12 'Vi + -^ 3 ^ 72 ' and I^rjii + contain 

only the noise from laser ^ 1 and have the following forms 


l2'ViLV2'-,3 — — I 1 I 2 ') Pi , (29) 

-^3^1' + ?73;2' = {'L>2''D2 — h'h) Pi . (30) 

Since in the quasi-stationary case any pairs of these operators commute, i.e., 'Di'Dji—'Dji'Di k, 
0, from Eqs. (12^ and [3(111 it is easy to derive the following expression for the unequal-arm 
interferometric combination that eliminates pp. 

~ h'h] (h'Vl + ^2';3) ~ [(^3^3' ~ -^1-^2')] ihVl' + ^3;2')- (31) 

If, on the other hand, the time-dependence of the delays is such as to prevent the delay oper¬ 
ators to commute, the expression of the unequal-arm Michelson combination above no longer 
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cancels pi. In order to derive the new expression for the nneqnal-arm interferometer that 
acconnts for “flexing” and simnltaneonsly cancels the noise from the microwave freqnency 
signal referenced to the onboard lasers, let ns hrst consider the following two combinations 
of the one-way measnrements entering into the observable given in Eq. flST]) : 

[/i/2'(/3?7i' -|- P3-2') + {h'Vl + h2';3);22'] = [D2'D2D^Dy — lilyI 2 'I^] Pi , (32) 

[-^l'-^3(-^2'hl + h2';3) + ihvi' + h3;2');3'3] = [-D3-D3'-D2'-D2 — -fl-^l'-^2'-^3] Pi • (33) 

Using Eqnations (1521 155]) we can use the “delay technique” again to finally derive the 
following expression for the new unequal-arm Michelson combination Xi^^. This new TDI 
combination accounts for the flexing effect and cancels the noise of the microwave signal 
referenced to the laser frequency through the OFC technique 

X^FC _ — I^I^,l2'h] [h'hih'Vl + V2'-,3) + {hVl' + P3;2');3'3] 

— [D^DyD 2 'D 2 — Iilyh'Is] [Iih'ihVl' + V3;2') + {h'Vl + h2';3);22'] • (34) 

As usual, and are obtained by cyclic permutation of the spacecraft indices. 

This expression is readily shown to be laser-noise-free to hrst order of spacecraft separa¬ 
tion velocities Lj (it is “hex-free”) and it simultaneously removes the phase noise from the 
microwave signal (generated by the OFC technique) used in the heterodyne measurements. 


D. The Sagnac combinations 


In the case of the Sagnac variables (a,/5,7,C) [l9| light originating from a spacecraft is 
simultaneously sent around the array on clockwise and counter-clockwise loops, and the two 
returning beams are then recombined. If the array is rotating and hexing, the two beams 
experience a diherent delay due to the rotation of the array (the Sagnac ehect) as well as 
the time-dependence of the armlengths and the resulting Doppler ehects. 

Since the previously derived Sagnac TDI combinations no longer simultaneously cancel 
the laser and the microwave signal noises, in order to derive the new TDI expressions let 
us hrst write down the six terms entering, for instance, into the a combination suitable 
for a stationary array (see equation (60) in 19|]) in the attempt to explicitly identify the 
expression of the new TDI combination, . This combination must satisfy the property of 
removing the same laser huctuations ahecting two beams that have been made to propagated 
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clockwise and counter-clockwise around the array. The expression for a (equation (60) of 
[i^ ) contains six terms, i.e. r]i',V 2 'r]^i,T>ii'D 2 iri 2 ',rii,V^r] 2 ,T>iD^ri^, which now assume the 
following forms 


hi' 

= T^2'P3 - h'Pi , 

(36) 

V2'r)3' 

= V 2 ''Dirp 2 — I^iV2'P3 , 

(36) 

'Dii'D 2 'T] 2 ' 

= T)iiT)2'^3'Pi ~ l 2 'T^vT^ 2 'P 2 ) 

(37) 

hi 

= T>3P2 - hpi , 

(38) 

T>3h2 

= — l2T>3,P2 , 

(39) 


= V 3 VIV 2 PI — J 3 D 3 D 1 P 3 . 

(40) 


By simple inspection of the above expressions it is easy to see that the following two linear 
combinations of the above six measurements only contain the laser noise pi 


= l2'l3'Vl' T -^2'^^2'h3' + 'DiiT) 2'V2' — [Dl''D2''Dy — , 

aj. = hhVi + hT^3V2 + — hhlslPi ■ 


(41) 

(42) 


Since the delay operators do not commute, we conclude that the straight difference at — a 


does not cancel the laser noise. However, if we now apply the “delay technique” 
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19| 


to the combinations given in Eqs. flTTl 1421) we hnally get the expression for the Sagnac 
combination 


a 


OFC — 


= — -^l-^2-^3]«t “ [Dv'D2i'D^i — Ji//2'/3/]a^ , 


(43) 


Although the combination still shows the presence of a residual laser noise, it can be 

shown to be small as it involves the difference of the clockwise and counter-clockwise rates of 
change of the propagation delays on the same circuit. For LISA, the remaining laser phase 
noises in i = 1, 2, 3, are several orders of magnitude below the secondary noises. 

In order to derive the new TDI expression for the fully symmetric Sagnac combination, 
C, we remind the reader that the rotation of the array breaks the symmetry and therefore 
its uniqueness. However, there still exist three generalized TDI laser-noise-free data combi¬ 
nations that have properties very similar to (C, and which can be used for the same scientihc 
purposes 


32j. These combinations, which we call )C 3 ^^)) can be derived by 


applying again our time-delay operator approach. 
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In order to proceed with the derivation of we should hrst remind ourselves that this 

TDI combination contains all six rji, rjii , i = 1, 2, 3, each of which is delayed only once 31 |. 
By simple inspection of the TDI combination (^i given in Eq. (67) of reference [l9j, we will 
start by considering the following six terms entering in it 


h3.3 = , 

h3',3 = T>3^^1'P2 — , 

Vi,i' = T^i'T^zPi - h'Di'Pi , 
r]r,i = ViV2'Pz - h'T^iPi , 

772,2' = T>2'T’iP3 — I2T^2'P2 , 
h2',2' = T^2'T^3'Pi ~ l 2 'T^ 2 'P 2 ■ 

( 44 ) 

If the array is rigidly rotating (i.e. its arm lengths are constant), it is easy to see that the 
following combinations contain only the noise from laser 1 

-^ 3 ^ 1 , 1 ' — hV3',3 + h'V3,3 = [l3'T^3T^2 — hh'I^i']pi , ( 45 ) 

h'Pl',! — h>V2,2' + 727 / 2 ', 2 ' = [I2'I^2''I^3' “ h'h'T^ilPi , ( 46 ) 


as we have used the commutativity property of the delay operators in order to cancel the p 2 
and ps terms. Since both sides of the two equations above contain only the pi noise, Ci^'" 
is found by the following expression: 


(^OFC _ _ jyJ 2 , Di \ (13771,1/ — 73773/,3 + 73/773,3) 


— [h'D3D2 — I 1 I 3 D 11 ] ( 72 / 771/,1 — 72 / 772 , 2 / + 727/2/, 2 /) , 


(47) 


which is a generalization of Eq. (67) given in 1^. If the delay-times also change in time, the 


perfect cancellation of the laser noises is no longer achieved in the i ( 2 ^^^ com¬ 

binations derived above. However, following the considerations made for the corresponding 
second-generation TDI combinations derived in ISj, it can be shown that for a mission 
like LISA the magnitude of the residual laser noises in these combinations are signihcantly 
smaller than the secondary system noises, making their effects entirely negligible. 
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E. The Monitor, Beacon, and Relay Combinations 


The expressions of the Relay Beacon and Monitor combinations 

that account for the rotation, flexing of the array, and implementation of the OFC technique, 
can also be derived by applying the time-delay iterative procedure highlighted in the previous 
subsections. However, as also noted in 1^, in order to reduce the number of terms in the 
monitor combinations we will rely on the use of the inverse delay operator D~^ , D~,^ , i = 
1,2,3. 


1. The relay 


In a relay combination, , one spacecraft is capable of only receiving a laser beam 

along one arm and transmitting along the other. If, for instance, we take spacecraft ^ 1 
to be the relay__spacecraft, the resulting TDI combination may contain the following four 


measurements 


18| : rjii , 7^2 , 7^2' , Vs'- If consider the expressions given in Eqs. (44- 


45) derived in 18| for the relay combination without OFC, it is relatively easy to obtain from 
them the following equivalent expressions valid when the OFC technique is implemented 


h'hh'V^' + hh'Vl',3' + l2VS'-,2'S' + V2;l'2'3' — h'l2'h'V2 — [T>3/'D2 ''Dii — Iyl2iIy]ViP3 ( 48 ) 

= T’l [T’l'T’3'’f^2' “ -fl'-f2'-f3']P3 (49) 

By applying the operator Vi[DiiV 3 iV 2 ' — I 11 I 21 I 31 ] to the left-hand-side of Eq. fHS]) . the oper¬ 
ator [P 3 /'D 2 '^^i' ~ Ivl 2 'l 3 ']T^i to the left-hand-side of Eq. fl49]) . and then take the difference 
of the resulting two expressions, we get the following expression for the second-generation 
relay combination when OFC is implemented 

= 'Di\DiiT>3iT>2i — Iiil2'l3'] {Ivl2l3'V2' + hh'Vl',3' + hVS'fl'S' + h2;l'2'3' “ h'l2'l3'V2} 
— ['D3i'D2''Dii — Iyl2'l3']'^l {h'V2'-,l'l + h'h'Vs',! + hl';3'l'l} (50) 

with 1 / 2 ^^ , obtained by cycling the spacecraft indeces. It can readily be verihed that 

the laser noise remaining in this combination vanishes to hrst order in the inter-spacecraft 
velocities. 
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2. The beacon 


In the beacon combination, say , spacecraft ^ 1 transmits (only) to the other two 
while the other two exchange laser light as usual. In this operating mode only the measure¬ 
ments t ]2 , r] 2 ' , r ]3 , r] 3 i will be available for synthesizing the combination . To identify 
its expression we proceed by modifying Eqs. (49-50) given in [l^ in such a way to exactly 
cancel the p 2 , Ps laser noises while retaining pi. After some long but straightforward algebra 
we get the following two expressions that contain only pi 

l2hvy,3' + hV2-,l'3' + h3;ll'3' ~ hl3'V3,3' = ’P3'[^l'^^l ~ -f2-f3']^^2Pl 5 (51) 

l2>h'V2,2 + h>V3'-,12 + V2'-,l'12 — hh'V2',2 = T^2[Dl'Dii — l2l3i]T>3'Pl . (52) 

By now applying our time-delay iterative procedure we obtain the hnal expression for P^^^ 

Pl^^ = T> 2 [T^iT^ 1 ' — hh']T^3'[hhV3',3' + hv2;l'3' + h3;ll'3' “ -^2-f3'h3,3'] 

— I)3i[DiiVi — l2l3']T^2[l2'I3'V2,2 + I2'P3'\12 + h2';l'12 “ hh'PT ,2\ (53) 


3. The monitor 


In the monitor combination, for instance, spacecraft ^ 1 receives (only) from the 

other two while the other two exchange laser light as usual; therefore relies only on 

the measurements rji , pi/ , r ]2 , r]y. If we now consider the expressions given in Eqs. (54-57) 
of ISj for the monitor combination without OFC, it is relatively easy to derive from them 
the following equivalent expressions valid when the OFC technique is implemented 


hl;ll' ~ hh'Vl = — l2l3']T^3P2 — — I2l3']pi , (54) 

-^3'h2,3 + h3';13 = ^^3 [^1^1' ~ -^2-f3']P2 ; (55) 

~ hh'Vl' = ~ -f2-f3']^2'P3 ~ ~ -^2-^3']^! ) (56) 

hvy,2' + V2-,l'2' = T^2'[T^1>T^1 — l2l3']P3 ■ (57) 


The above expressions can be hrst combined in pairs to remove the p 2 , P 3 noises in the 
following way 

[Dl'Vl — l2h']'^3[h'V2,3 + h3';13] ~ ~ -^2-^3'] [hi; 11' “ -f2-f3'hl] 

= -fl^^3[^l^l' ~ -f2-f3'][^^l'^l ~ -^2-^3']^! ) (58) 
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[DiDy — l2h']T^2'[hV3',2' + ^2;1'2'] “ X>2'“ -^2-^3'] [^l';l'l “ -^2-^3''^!'] 

= IiiV2'[DiiVi — l2ly][DiVii — l2l3']Pl ■ (59) 


To find the final expression for we could apply again the iterative procedure. However, 


to get an expression that has a smaller number of terms 
operators, By simply inspecting the right-hand-sides of Eqs. 

to derive the following expression for 


18[ |. we will use the inverse delay 

, it is easy 


E 


OFC 

1 


— lyV^ ^[DiiVi — l2h']T^3[h'V2,3 + h3';13] “ h'[DiDy — /2/3'][hl;ll' ~ -^2f^3'hl] 

-^1^2'^ [^1^1' “ -^2.f3']^^2'[.f2h3',2' + h2;l'2'] + h[T^l'T^l “ f^2f^3'] [hl'P'l ~ l2h'Vl']{Q^) 


which can easily be shown to be laser noise-free to hrst order in the systematic relative 
velocities of the spacecraft. 


IV. SUMMARY AND CONCLUSIONS 

We have derived second-generation time-delay interferometric combinations valid when 
the microwave signal used for heterodyning the phase measurements is generated by an 
onboard optical-frequency comb subsystem rather than a USO. This provides a microwave 
signal that is coherent to the frequency of the onboard stabilized laser, and results in a sig- 
nihcant simplihcation of the onboard interferometry system. This is because (i) generation 
of modulated beams and additional heterodyne measurements involving sidebands are no 
longer needed, and (ii) the entire onboard USO subsystem can be replaced with the mi¬ 
crowave signal referenced to the onboard laser. The corresponding hardware simplihcation 
results in a considerably reduced system complexity and probability of subsystem failure. 
Recent progress in the realization of a space-qualihed OFC indicates that such a capabil^ 
will be available before the planned hight of the ESA gravitational wave mission eLISA |6|. 
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